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Abstract. We establish a close quantitative analogy between the excitation and ionization process of highly 
excited one electron Rydberg states under microwave driving and charge transport across disordered ID 
lattices. Our results open a new arena for Anderson localization - a disorder induced effect - in a large 
class of perfectly deterministic, decaying atomic systems. 

PACS. 32.80.Rm Multiphoton ionization and excitation to highly excited states (e.g., Rydberg states) - 
05.45.Mt Semiclassical chaos / quantum chaos - 72.15.Rn Localization effects (Anderson or weak local- 
ization) 
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1 Introduction 

Probability transport in disordered media is an equally 
exciting and active field at the very heart of statistical 
physics. Many particle dynamics in a gas, turbulent hy- 
drodynamic flow, the traffic flow across canonically over- 
crowded european highways, and, last but not least, the 
flow of money across stock and option markets are de- 
scribed by statistical means. Also in the microscopic realm 
we are often confronted with similar situations where it is 
impossible to track each single detail of the system under 
study - either because we ignore the precise form of the 
potentials which generate the dynamics, or because the 
number of individual constituents of the system is simply 
too large, or a combination of both. Well-known exam- 
ples thereof are the charge transmission through meso- 
scopic wires, the scattering of slow neutrons off heavy nu- 
clei, or the ion transport across cellular membranes. On 
this microscopic level, however, when noise is sufficiently 
weak, quantum interference and tunneling can dramati- 
cally affect classical probability transport, and weak 
and strong (vulgo Anderson) localization the Mott- 
Hubbard metal-insulator transition (see ZoUer's contribu- 
tion to this issue), Ericson fluctuations j^, chaos assisted 
tunneling Q, or (universal) conductance fluctuations Q 
are just the most prominent ones of the many surpris- 
ing phenomena which are born out by coherent complex 
transport on the microscopic scale. 

On a first glance, such signatures of complex dynam- 
ics (brought about by random potentials and/or many- 
particle interactions) appear precisely what we do not ex- 
pect when dealing with well isolated atomic or quantum 
optical systems. In the quest for an almost perfect control 



of matter, we seek optimal control of the potentials which 
determine the dynamics (see Gerber's contribution), in 
time and in space. Hence, how can such low-dimensional 
systems mimic the complex dynamics of disordered sys- 
tems? The answer is - precisely - Hamiltonian chaos, in- 
duced by nonlinear and perfectly controlled coupling of 
the few (at least two are needed Q) degrees of freedom 
available. Given a sufficiently high density of states (in 
a quantum system with a discrete or quasi-discrete spec- 
trum), such coupling can destroy the symmetries, and this 
is the good quantum numbers of the unperturbed dynam- 
ics, giving rise to an extremely complex spectrum charac- 
terized by an abundance of anticrossings of strongly vari- 
able size (see also Saenz' contribution) 0. If we now pre- 
pare a wave function of arbitrary shape at an arbitrary lo- 
cation in phase space, its time evolution will indeed reflect 
the complex spectral structure, and, consequently, feature 
the characteristic properties of complex/disordered sys- 
tems. This is the essential consequence of the random ma- 
trix conjecture which states that the spectral prop- 
erties of low dimensional quantum systems with underly- 
ing chaotic classical dynamics exhibit the same (universal) 
statistical features as complex quantum systems described 
by Random Matrix Theory (which took its origin in the 
attempt |^ to give a robust description of compound nu- 
clear reactions - despite the little knowledge about the 
detailed nature of the potentials which generate the ex- 
perimentally observed cross sections). 

Which are the good experimental observables to mon- 
itor those features? In the usual setting of, e.g., charge 
transmission across a disordered lattice, transport occurs 
in conflguration space, and the experimentalist measures 
a transmission probability. Such scenario can nowadays 
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be mimicked in quantum optical table top realizations of 
"ideal" lattices, where a suitable laser configuration es- 
tablishes a periodic (ID or 2D or 3D) potential, with 
adjustable lattice constant and modulation depth. Cold 
atoms or ions initially prepared at one specific location can 
then be monitored as they move across the sample | p^ , 
and it is equally possible to enforce the transition from lo- 
calized to delocalized (in configuration space) eigenstates 
on the lattice On the other hand, we may ask with 
equal right whether certain phenomena familiar from the 
theory of complex and/or disordered systems can be im- 
ported to electronic dynamics on the scale of a single, 
simple atom, with few degrees of freedom, and no a pri- 
ori resemblence with familiar transport problems. In other 
words, do disorder induced phenomena have any general 
and robust bearing for the dynamics of simple atomic 
systems exposed to strong perturbations? The answer is 
positive, and can be elaborated for various of the above- 
mentioned coherent transport phenomena. For reasons of 
time and space, we shall specialize here to Anderson local- 
ization in driven atomic systems, and, more specifically, in 
atomic Rydberg states under microwave driving. 



2 Anderson localization on the atomic scale 

Anderson localization is a disorder induced effect, which 
was predicted |^ by Anderson to occur in the transmis- 
sion of a charge across a one dimensional, disordered lat- 
tice. As depicted in Fig. |l|, the problem cooks down to 
the transmission of a particle across a one dimensional 
random potential, at a given injection energy. The num- 
ber of paths which guide the particle from the left to the 
right of the sample is exponentially large, since at each 
hump of the potential the particle can be either reflected 
or transmitted (with randomly distributed transmission 
and reflection coefficients), thus multiplying the available 
paths by a factor of two at each hump. The transmission 
amplitudes associated with these different paths acquire 
randomly distributed phases as they migrate through the 
sample, and interfere destructively (if the phases get ho- 
mogeneously distributed over the unit circle) on exit at 
the right edge of the sample. Hence, a disordered potential 
leads to the suppression of transmission by quantum inter- 
ference, and a more quantitative analysis shows that the 
eigenfunctions of the particle are exponentially localized 
over the lattice, with a characteristic localization length ^. 
The final figure of merit which determines the measured 
conductance across the sample is £,/L, with L the length of 
the sample. ^/L determines the population of the last site 
right at the edge of the sample, and hence the probability 
flux which can escape from the sample, via the coupling 
matrix element connecting the last site to the leads - at- 
tached to the sample to probe the conductance. 

We now import this scenario to the atomic realm, more 
specifically to the excitation and ionization dynamics of 
atomic Rydberg states under microwave driving. As men- 
tioned above, for features of complex/disordered trans- 
port to become manifest on the atomic scale, we need 
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Fig. 1. The Anderson scenario: A particle (dashed horizontal 
line) is transmitted across a disordered potential in ID config- 
uration space. At each hump, the particle can be reflected or 
transmitted, with random reflection/transmission coeflicients. 
The atomic ionization problem under microwave driving at fre- 
quency u is strictly analogous if we replace reflection and trans- 
mission by emission and absorption of a photon into/from the 
driving fleld, and substitute configuration space by the energy 
axis. The "atomic sample length" L is then the ionization po- 
tential expressed in integer multiples of the photon energy u, 
and the "sample edge" is defined by the continuum threshold. 
The atomic decay to the continuum (indicated by the right- 
most arrow across the ionization threshold) is mediated by the 
one photon transition matrix element from the highest bound 
Rydberg state - connected to the initial state \(j)o) (indicated 
by the vertical arrow) by a sequence of near resonant (though 
slightly detuned - small verticle ticks along the horizontal en- 
ergy axis indicate the near resonantly coupled bound states, 
dashed verticle lines highlight the lattice period) one photon 
transitions. The population of this last bound state depends 
exponentially on L/^, with ^ the characteristic decay length of 
the localized wavefunction. 



a high density of states actively involved in the dynam- 
ics. This is guaranteed in the Rydberg regime under mi- 
crowave driving, since the typical Rydberg energy split- 
ting lies in the microwave frequencies, and, consequently, 
many Rydberg states will be efficiently coupled by the 
drive (through subsequent, near- resonant one photon ab- 
sorptions). The atomic transport process which we want 
to parallel with the Anderson scenario is the transport of 
electronic population initially prepared in a well defined 
atomic initial state \4>o) = |no, Iq, mo) (where spherical 
quantum numbers are used, as will be done throughout 
the sequal of this paper) towards the atomic continuum 
- easily measured as the ionization yield Pjon for given 
fleld amplitude F, field frequency w, and interaction time 
t. To make the analogy complete, we have to remind our- 
selves of the slow variation of the level splitting with en- 
ergy in the Rydberg domain. As a consequence of this 
unharmonicity of the spectrum, the subsequent one pho- 
ton transition matrix elements which establish the "hop- 
ping matrix elements" between the different "sites" on the 
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energy axis mimic a random series , due to the imper- 
fect matching (tantamount to quasi random detuning) of 
the driving frequency with the exact transition frequency 
(much as the simplest random number generators built 
on a modulo operation Q ) . Consequently, if we replace 
"transmission/reflection" in the above transmission prob- 
lem a la Anderson by "absorption/emission" in the atomic 
excitation process, ID configuration space by the energy 
axis, and the attached leads by the atomic continuum, we 
encounter precisely the same scenario: An exponentially 
large number of paths visiting different sites (near res- 
onantly coupled Rydberg states in the atomic problem) 
accumulate quasi random phases on exit from the sam- 
ple, i.e., at the ionization threshold. If one photon tran- 
sitions from the last near resonantly coupled bound state 
to the continuum define the dominant ionization chan- 
nel, the atomic ionization yield will be proportional to 
this state's population - determined by the localization 
parameter ^/L already familiar from the Anderson sce- 
nario. Again, the sample length L is given by the number 
of sites which span the lattice, equivalent to the ioniza- 
tion potential of the atomic initial state in multiples of 
the driving field photon energy. Provided ^/L ^ 1, the 
Anderson model predicts exponential suppression of the 
ionization yield due to quantum interference. Despite the 
strong, near resonant driving, the atom will cease absorb- 
ing energy from the field! 

In fact, there already are abundant experimental data 
which indicate that the above mechanism 
is at work in driven Rydberg states of atomic hydrogen 
and of alkali atoms. The central experimental result is 
an increase of the "scaled ionization threshold amplitude" 
^^(10%) _ ^(io%)j^4 ^\^]^ increasing principal quantum 
number rtg of the initial atomic state. In other words, 
the field amplitude F^'^''"' required to ionize 10% of the 
atoms, at given field frequency oj and interaction time t, 
measured in units of the average Coulomb field tiq"^ 
experienced by the electron on its unperturbed initial Ke- 
pler orbit, increases as we decrease the ionization potential 
of the atomic initial state! This is in dramatic contradic- 
tion with the result of a classical treatment, which pre- 
dicts systematically smaller ionization thresholds (due to 
classically chaotic phase space transport) than observed 
in the experiment. Consequently, there exists a range of 
field amplitudes which induce efficient classical ionization, 
whilst the experimentally observed yield is close to zero. 
This quantum suppression of classically chaotic ionization 
is interpreted as a signature of Anderson localization, and 
traditionally termed "dynamical localization" , such as to 
identify dynamical chaos as the cause of the quasi-rando- 
mization of the hopping matrix elements in the Anderson 
picture. However, this experimental finding is only con- 
sistent with the hypothesis of Anderson localization, it 
is not a clear proof - simply since one may imagine dif- 
ferent mechanisms which inhibit the ionization process, 
such as semiclassical stabilization effects caused by (par- 
tial) barriers in phase space ||l^,|l^,|o|. Furthermore, most 
experiments have been performed under slightly different 
experimental conditions: not only have different atomic 



species been used, but also different atomic initial states 
were exposed to microwave fields of different frequencies u 
and of variable duration t. All available experimental data 
exhibit large quantitative differences between the ioniza- 
tion threshold of atomic hydrogen on the one hand and of 
nonhydrogenic initial states of alkali atoms on the other, 
with the alkali thresholds down by a factor five to ten as 
compared to the hydrogen thresholds [|9|. The original 
theory of dynamical localization in driven atomic systems 

, based on a very crude model of the actual bound state 
atomic dynamics, is not suited to explain these quantita- 
tive differences nor to definitively exclude processes dis- 
tinct from Anderson localization which might cause the 
observed increase of . 

What we therefore need is an accurate theoretical treat- 
ment of microwave driven one electron Rydberg states of 
atomic hydrogen and of alkali atoms. Such treatment is 
not only required to explain the experimental findings so 
far available in a consistent and unified way, but also to 
guide future experiments which seek to test quantitative 
predictions which follow from Anderson localization the- 
ory for the atomic ionization process. 



3 An accurate treatment of one electron 
Rydberg states under electromagnetic driving 

Our specific problem is described by the 27r/w-periodic 
Hamiltonian 



H{t) = 



V;tom(r) + F • r coscjt, r > 0, (1) 



where Vatom('') denotes the atomic potential seen by the 
valence electron, which is taken care of by a R-matrix 
approach |Q |^ The periodicity of H suggests to explore 
Floquet theory |24[^, and we therefore end up solving 
the Floquet eigenvalue problem 



with the Floquet Hamiltonian 

n = H -idt, 



(2) 



(3) 



the spectrum of which is invariant under translations by 
to. Knowledge of the \ej) and ej within an energy range of 
width Lo is therefore sufficient for a complete description 
of the dynamics. 

After a further Fourier transform of the 27r/cL;-periodic 
\e.j ) , amended by complex dilation of the Hamiltonian [ p5| 
p6| , the eigenvalue problem is represented in a real Stur- 
mian basis, what converts (^ into a block-tridiagonal. 
cornplex symmetric, sparse banded eigenvalue problem p2[ 
Note that the strong selection rules induced by the 
Sturmian basis are absolutely crucial for the numerical 



^ The R-matrix treatment requires input of the quantum de- 
fects of the non-hydrogenic angular momentum states of the 
unperturbed alkali atoms, known from high precision spec- 
troscopy experiments [E3|. 
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treatment of our problem, since they allow for a consider- 
able reduction of the required memory. In the parameter 
range typically employed in the laboratory, typical dimen- 
sions are 10^ x 10^, what nontheless requires a very large 
parallel supercomputer. Indeed, all the results presented 
hereafter were obtained on the CRAY T3E of the High 
Performance Computing Center RZG of the Max-Planck- 
Society at Garching and, once this machine became too 
small, on the Hitachi SR8000-F1 of the Bavarian Academy 
of Sciences in Munich. Though, the availability of such a 
large machine is not enough. In addition, an efficient paral- 
lel implementation of the Lanczos diagonalization routine 
is needed, which is by no means a trivial requirement, due 
to considerable communication between different (and not 
only adjacent) processors of the parallel machine. 

Once the theoretical and numerical machinery briefly 
sketched above is assembled, we can start our "numeri- 
cal experiment" , which closely mimics the reality in the 
laboratory. The "counts" which the numerical experimen- 
talist has to collect much as his colleague in the real lab 
are the poles of the resolvent operator in the lower half 
of the complex plane [^ , i.e. the complex eigenvalues 
Ej = Ej — of the above eigenvalue problem. From 

these we can immediately extract the ionization yield psf 

Pio„ = l-5]|(0o|e)|'exp(-r,t), (4) 

e 

for the specific choice of |0o): and F. Note that ap- 
prox. 4000 . . . 10000 poles are typically collected here. This 
is simply due to the fact that the decomposition of the 
atomic initial state |0o) in the dressed state basis is ex- 
tremely broad in our present case [p9| , what is a direct 
consequence of the efficient destruction of good quantum 
numbers by the driving microwave field, and, hence, of 
quantum chaos. 



4 A unified view on Anderson localization in 
driven atoms 

With the above, we can now perform a direct compari- 
son of the ionization thresholds of atomic Rydberg states 
of hydrogen and of those of alkali Rydberg states, under 
precisely equivalent conditions. For the sake of comparison 
to the arguably largest experimental data set (produced 
in the Stony Brook group Q), we specifically choose 
w/27r = 36 GHz, i = 327 x 2tt/u, and = 28... 80, 
^0 — ™o = 0, where mo is conserved under linearly po- 
larized driving. Precisely as in the laboratory, we scan F 
from low to high values, monitor the dependence of Pion on 
F, and extract Fq^"^"-*, for increasing values of uq. Fig. || 
shows a comparison of our numerical results for lithium 
£o = states to laboratory results on atomic hydrogen. 
Surprisingly, and in perfect contradiction to all published 
experimental evidence, the nonhydrogenic (the = 
state carries the largest quantum defect 5i„=o = 0.399468) 
lithium Rydberg states exhibit essentially the same ion- 
ization thresholds as atomic hydrogen in regime (I), i.e. 




Fig. 2. Comparison of the ionization threshold of atomic hy- 
drogen (obtained in laboratory experiments |l4|^ ; light sym- 
bols) and of lithium (^o = 0) Rydberg states (from our nu- 
merical experiment; dark symbols) under microwave driving 
of duration t = 327 x 2tt/u> at uj/2tt — 36 GHz, in the range 
no = 28 ... 80 of principal quantum numbers (from left to right, 
since loq = ujuq). Since the laboratory values of the physical 
parameters are identical, the present plot in scaled units does 
not imply any a priori assumptions on alkali scaling laws! 

for "scaled frequencies" luq = ujuq > 0.8. Only for de- 
creasing values of no leading to smaller values of loq do 
we observe an increasing discrepancy between the hydro- 
genic and alkali thresholds. It turns our that this can be 
traced back to the local density of states in the Rydberg 
series of the different atomic species : For to larger than 
the local hydrogenic level spacing - which scales like , 
hence loq > \ - the atomic excitation process is essen- 
tially unaffected by the spectral substructure in the al- 
kalis. The atomic system offers a ladder of near resonant 
one photon transitions connecting |(/)o) to the continuum 
(see Fig. 0), and it is irrelevant for the transport process 
whether this occurs under the participation or in the ab- 
sence of non-hydrogenic states. However, once the driving 
frequency is smaller than the hydrogenic level spacing but 
still larger than the splitting between the non-hydrogenic 
alkali initial state and the hydrogenic manifold, the same 
sequence of near resonant one photon transitions can still 
be established in the alkali atom, whilst in hydrogen only 
sequences of higher order transitions are left to mediate 
the ionization process. Hence, the realm of the Ander- 
son scenario outlined in the introduction extends over 
a wider range of principal quantum numbers in the al- 
kali species than in atomic hydrogen, and this is indeed 
manifest in Fig. |^: In regime (II), with 0.4 < wo ^ 0.8, 
the lithium data continue to exhibit decreasing thresholds 
with decreasing woj ^ an indicator of Anderson localiza- 
tion, whilst the hydrogenic thesholds rapidly increase as 
the binding potential is increased (and the local density 
of states is decreased, as compared to the photon energy). 
Finally, in regime (HI), also the local energy splitting in 
the alkali atom is too large to be driven near resonantly 
by a single photon. The Anderson scenario breaks down, 
and the thresholds saturate in the plot of Fig. ^. 

The above results hold, at least, a two-fold message: 
First, there is the prediction that comparable ionization 
thresholds for alkali and hydrogen Rydberg states will be 
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observed once the driving frequency is larger than the 
spacing between adjacent hydrogenic manifolds. This quan- 
titative prediciton can be checked immediately in state 
of the art laboratory experiments [p^ . Second, when re- 
analyzing all available experimental data on the microwave 
ionization of non- hydrogenic alkali Rydberg states by scal- 
ing the driving field amplitude and frequency as done in 
Fig. H, we find |^ that all these data were obtained in 
regimes (II) and (III), what explains the apparently dra- 
maticly enhanced ionization yields of non- hydrogenic al- 
kali Rydberg states as compared to hydrogen |^,^. Note 
that this discrepancy remained a puzzle over more than 
one decade |Q, since, so far, precisely identical experi- 
mental conditions for different atomic species were never 
established in the lab, and the appropriate scaling of F 
and uj was controversial due to the badly defined classical 
analog of the alkali Rydberg dynamics As a 

matter of fact, our comparison of hydrogen and lithium 
data in terms of scaled units in Fig. g do not imply any a 
priori assumptions, since they were obtained for the same 
laboratory values of all physical parameters, but the fact 
that alkali and hydrogen thresholds coincide in regime (I) 
proves that, for sufficiently high driving frequencies, the 
hydrogenic sca ling holds even for alkalis! Even more, we do 
know by now |3l| that also non-hydrogenic initial states 
of sodium and rubidium exhibit the same, hydrogenic, 
threshold in regime (I), what strongly suggests that this 
is a universal threshold for one-electron Rydberg states. 
Only this universality makes the increase of p^'^^'^"^ with 
ujo or rig a sufficient condition for the Anderson scenario 
to prevail, since only the universality shows that the only 
relevant ingredient is the local density of states compared 
to the photon energy. 

There is another consequence of Anderson localization 
theory which can be imported to our atomic system: The 
localization length ^ is really well-defined a quantity only 
in the limit of infinite sample lengthy and will fluctuate 
around this limiting value for finite L |32| . The theoretical 
modelling of the conductance across an Anderson localized 
ID wire of finite length L predicts a normal distribution of 
the logarithm of the properly normalized conductance g, 
if sampled over different realizations of the random lattice 
potential (see Fig. at fixed ^/L |32 . Since the (meso- 
scopic) conductance is given by nothing but a sum over 
transition matrix elements (from left to right of the sam- 
ple) , we can come up with an analogous definition of 
the "atomic conductance" in our ionization problem pm , 
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(5) 



since the decay rates Fj can be understood as transition 
matrix elements of a suitably defined Floquet scattering 
problem (with A the average level spacing). Now, if 
Anderson is at work in the atomic problem, the atomic 
conductance has to exhibit the same statistical properties 
as required by the theory of disordered transport - and, 
indeed, it does! Fig. ^ shows the distribution of the atomic 
conductance of a ID model atom |Q (the electron is con- 
fined to the configuration space axis parallel to the polar- 




Fig. 3. Distribution (histograms) of the atomic conductance 
g of a ID Rydberg atom sampled over 500 differ- 

ent realizations of the localization parameter ^/L — 0.2, in 
the frequency range uo — 2.0 . . . 2.5, for no = 40 (left) and 
no — 100 (right). The lognormal fit (smooth curve) is excellent 
for no = 100, in perfect quantitative agreement with Ander- 
son localization theory |32]. The deviations from lognormal 
behaviour for no — 40 reflect a finite size effect. 



ization of the driving field), for fixed ^/L = 0.2 (which is 
given in terms of no, i^, w, according to the original theory 
of dynamical localization |^), and for two different val- 
ues uq = 40, 100 of the principal quantum number The 
lognormal fit to the data is excellent for uq = 100, what 
represents another, independent and quantitative indica- 
tor of Anderson localization as the dominant mechanism 
which determines the energy exchange between the atom 
and the field. We also observe in Fig. ^ that the fit is 
significantly worse for no — 40, what, however, does not 
contradict our preceding statement: At too low quantum 
numbers, the atomic sample size L, i.e. the number of near 
resonantly coupled Rydberg states between |(/)o} and the 
continuum threshold is too small (L ~ 10) to allow for 
a smooth exponential localization of the electronic wave- 
function over the energy axis (^ = 0.2L = 2). The figure 
therefore simply highlights a finite size effect. 

Finally, Fig. ^ shows the physical imprint of the lognor- 
mal distribution of Fig. |^, the dependence of the ionization 
yield Pion on the scaled frequency ujq, at ^/L = 0.2 (be- 
low threshold in Fig. H) : The yield is typically very small, 
close to zero, but exhibits strong, erratic enhancements 
at particular values of loq. Correspondingly, the atomic 
conductance g fluctuates on a logarithmic scale, in the 
right column of the flgure. This is nothing but the imme- 
diate signature of the highly sensitive interference of the 
many paths defined by subsequent absorption and emis- 
sion events mediating the transition from the initial state 
to the atomic continuum (see Fig. |l]), and a specific sig- 
nature of Anderson localization. The reader should con- 
template this strongly fluctuating signal, since what fluc- 
tuates is a quantity obtained from a weighted average over 
the entire Floquet spectrum, and not just one single rate 
r ' 



^ For each plot approx. 500 complex-valued spectra have 
been sampled, making such statistics even for the ID model 
atom a rather expensive enterprise. By now, however, we have 
first evidence from 3D calculations that the same result can be 
expected for the real atom |29(]. 
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Fig. 4. Ionization yield Pion (left) and atomic conductance 
g (right) of a ID Rydberg atom p^,^, as a function of the 
scaled frequency uio, in the Anderson/dynamically localized 
regime at ^/L = 0.2, no = 100 (below threshold in Fig. ^). 

5 Conclusion 

We provided two independent, quantitative proofs for An- 
derson localization as the dominant mechanism which gov- 
erns the excitation and ionization process of strongly driven 
one electron Rydberg systems. Given the universal ion- 
ization threshold which we observed for atomic hydrogen 
and lithium Rydberg states in the high frequency param- 
eter range, it appears legitimate to speculate that this 
disorder-induced quantum interference effect can be gen- 
eralized for an even larger class of driven atomic or molec- 
ular system, if only the simple Anderson scenario can be 
etablished on the energy axis, as depicted in Fig. [I]. 
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